APPROXIMATION FOR THE DISTRIBUTION OF EXTREMES 
OF ONE DEPENDENT STATIONARY SEQUENCES OF 
RANDOM VARIABLES 

AMARIOAREI ALEXANDRU 



Abstract. In this paper we improve some existing results concerning the ap- 
proximation of the distribution of extremes of a 1-dcpcndcnt and stationary 
sequence of random variables. We enlarge the range of applicability and im- 
prove the approximation error. An application to the study of the distribution 
of scan statistics generated by Bernoulli trials is given. 



1. Introduction 

The starting point of this paper is a series of results of Haiman [1999], concerning 
the extreme of a 1-dependent and stationary sequence of random variables. Let 
(X n ) n >i be a sequence of strictly stationary 1-dependent random variables (for 
any t > 1 we have o~(X\, . . . , X t ) and o~(X t +2, ■ • ■ ) are independent) with marginal 
distribution function F(x) = P(Xi < x). Let x such that 

ini{u\F(u) > 0} < x < sup{u\F(u) < 1}. 

Define the sequences 

Pn = Pn{x) = P(min{Xi, X 2 , . . . ,X n } > x), n > 1, p = 1, (1.1) 
q n = q n {x)=F(m&x{X 1 ,X 2 ,...,X n }<x), n>l (1.2) 



and the series 



C(z) = C x (z) -1 + J] (-l)V-iz' (1.3) 

k=l 

In Haiman [1999], the author proved the following results: 

Theorem 1.1. For x such that < pi(x) < 0.025, C x (z) has a unique zero X(x), 
of order of multiplicity 1, inside the interval (1, 1 + 2pi), such that 

|A - (1 + pi - p 2 + p 3 - Pi + 2p? + ip\ - 5pipa)| <87p?. (1.4) 

Theorem 1.2. We have 

gi = l— pi, q2 = 1 - 2pi +p 2 , 93 = 1 - 3pi + 2p 2 +p\ —p3 
and for n > 3 if pi < 0.025, 

\q n X n -(l-p 2 + 2p 3 -3p 4 +pl + 6p 2 2 -6 Pl p 2 )\ < 561p?. (1.5) 
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These results were successfully applied in a series of applications: the distribu- 
tion of the maximum of the increments of the Wiener process (Haiman [1999]), 
extremes of Markov sequences (Haiman et al. [1995]), the distribution of scan sta- 
tistics, both in one dimensional (see Haiman [2000, 2007]) and two dimensional case 
(sec Haiman and Prcda [2002, 2006]). 

Following the same lines of proofs as in Haiman [1999], we improve Theorem 1.1 
and Theorem 1.2 by enlarging the range of applicability and providing sharper er- 
ror bounds. The main results are presented in Section 2. In Section 3, we present 
an application to the study of the distribution of one dimensional discrete scan 
statistics emphasizing the difference between the new and the old results. Proofs 
are presented in Section 4. 

2. Main results 

The following theorem gives a parametric form of the Theorem 1.1 improving both 
the range of pi(x), from 0.025 to 0.1, and the error coefficient: 

Theorem 2.1. For x such that < pi(x) < a < 0.1, C x {z) has an unique zero 
X(x), of order of multiplicity 1, inside an interval of the form (1, 1 + Ipi), such that 

\\-(l+Pi-p 2 +P3-P4 + 2p 2 1 + 3p 2 2 -5p lP2 )\ <K{a)pl (2.1) 

where I = 1(a) > t^(a), £2 (a) is the second root in magnitude of the equation 
at 3 — t + 1 = and K(a) is given by 

ll-3a 1 97/1 , Q „ A 2+3Za-a(2-Za)(l+ia) 2 
[l-a(l+la) 2 ] 2 

Using the properties of the 1-dependent sequence, an immediate consequence of 
Theorem 2.1 is the following: 

Corollary 2.2. Let A be defined as in Theorem 2.1, then 

\\ - (1 + Pl - p 2 + 2{p x - p 2 f)\ < (1 + aK{a))p\. (2.3) 

To get a better grasp of the bounds in Theorem 2.1 and Corollary 2.2, we present, 
for selected values of a, the values taken by the coefficients in Eq.(2.1) and Eq.(2.3): 



Q 


I 


K(a) 1 - 


1- aK(a) 


0.100 


1.5347 


38.6302 


4.8630 


0.050 


1.1893 


21.2853 


2.0642 


0.025 


1.0835 


17.5663 


1.4391 


0.010 


1.0313 


15.9265 


1.1592 



Table 1. Selected values for the error coefficients in Theorem 2.1 
and Corollary 2.2 



Notice that for the value considered in Haiman [1999], i.e. a = 0.025, our corre- 
sponding value for the error coefficient in Eq.(2.1) is almost five times smaller than 
in Eq.(1.4). The following result improves Theorem 1.2: 
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Theorem 2.3. Lets suppose that x is such that < pi(x) < a < 0.1 and define 
i] = 1 + la with I = 1(a) > t%(a) and £2(0:) the second root in magnitude of the 
equation at 3 — t + 1 = 0. // A = X(x) is the zero obtained in Theorem 2.1, then the 
following relation holds 

\qnX n -(l-P2+2p3-3p i + pl + 6p 2 2 -6p 1 p 2 )\<r(a)p 3 1 , (2.4) 

where T(a) = L(a) + E{a), K(a) is given by Eq.(2.2) and 

L{a) = 3K{a){l + a + 3a 2 )[l + a + 3a 2 + K(a)a 3 } + a 6 K 3 (a) 

+9a(4 + 3a + 3a 2 ) +55 (2.5) 

pM nil V 5 [j + iiz 2q ^ 4 t 1 + "fa ~ 2)] [1 + g + (1 - 3a)7? 2 ] 

W ' 2(l-a77 2 ) 4 [(l-a77 2 ) 2 -a77 2 (l + 77-2ar ? ) 2 ] 1 ' 

The next corollary is an immediate consequence of Theorem 2.3: 

Corollary 2.4. In the conditions of Theorem 2.3 we have 

\q n \ n -(l-p 2 )\ < (3 + aT(a))pl (2.7) 

The error coefScient in Eq.(2.4) is smaller in comparison with the corresponding 
one from Eq.(1.5), as the following table can show: 



a 


T(a) 


3 + «r(a) 


0.100 


480.696 


51.0696 


0.050 


180.532 


12.0266 


0.025 


145.202 


6.6300 


0.010 


131.438 


4.3143 



Table 2. Selected values for the error coefficients in Theorem 2.3 
and Corollary 2.4 



Combining the results obtained in Theorem 2.1 and Theorem 2.3 we get the fol- 
lowing approximation: 

Theorem 2.5. Let x such that qi(x) > 1 — a > 0.9. IfT(a) and K(a) are the 
same as in Theorem 2.3, then 



6(91 - q 2 ? + Ms - 3q 4 



(1 + qi - q 2 + 93 - 94 + 2qf + 3q 2 - 5 qi q 2 ) n 



<Ai(l-gi) 3 (2. 



with 



Ai = A 1 (a,n)=T(a)+nK(a). (2.9) 

In the same fashion combining the results from Corollary 2.2 and Corollary 2.4 we 
get 

Theorem 2.6. If x is such that q\[x) > 1 — a > 0.9, then 



2?1 - 92 



[1 + 91 -92 + 2(9i -92)T 



< Aa(l-ft) 



with 



A 2 = A 2 (a, n, 91) = 3 + T{a)(l - qi ) + n [1 + K(a)(l - 91) 
and where T(a) and K(a) are the same as in Theorem 2.3. 



(2.10) 
(2.11) 
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3. Application to the distribution of scan statistics 

In many applications the decision makers have to determine if a certain accumu- 
lation of events is normal or not, where by normal we mean that it can be ex- 
plained by an underlying probability model defined by a null hypothesis of ran- 
domness. One way of determining if such a cluster of observations is exception- 
ally or not is by using the scan statistics. For a comprehensive methodologi- 
cal treatment and a rich source of applications of this test statistic, one may 
study the books of Glaz, Naus and Wallenstein [2001] and more recently the one of 
Glaz, Pozdnyakov and Wallenstein [2009]. 

Let Y\,Y2, . . . , Y/v be a sequence of independent and identically distributed (i.i.d.) 
random variables and m. 1 < m < N, be a fixed positive integer. If we consider 
the random variables 

t+m-1 

Z t = ]T Yi, \<t<N -m + l (3.1) 

then the one dimensional discrete scan statistic is defined by 

S m (N) = max Z t . (3.2) 

l<t<N-m+l 

To get an intuitive meaning of the above definition, if Yi are integer valued random 
variables then we can interpret them as the number of observed events at the time 
i. The scan statistics is then viewed as the maximum number of events observed 
in any contiguous period of length m within the interval {1,2,..., N}. Since exact 
formulas for the distribution of S m exist only in a small number of situations (see 
for example Glaz, Naus and Wallenstein [2001], Chapter 13), various approxima- 
tion methods and bounds have been proposed. In what follows we will use the 
approximation method developed in Haiman [2007] but with the help of the results 
obtained in Section 2. The method is based on the important observation that, 
in the i.i.d. case, the discrete scan statistic can be expressed as an extreme of a 
1-dependent stationary sequence. It is easy to see that the random variables 

W k = max Z s , fc = l,2,... (3.3) 

(k-l)m+l<s<km+l 

form a 1-dependent stationary sequence and that the following relation holds, 

S m {Lm) = max (W u W 2 ,..., W L -i). (3.4) 

Observe that for general N, one can consider L = [N/m] and then apply the 
inequality 

P(S m [(L + l)m] <n) <¥(S rn (N) < n) <P(S m (Lm) <n). (3.5) 

When P(Wi > n) < 0.1, we can apply the results from Theorem 2.6 to the sequence 
Wi, W2, ■ ■ ■ , Wl-i to obtain the following approximation for the distribution of the 
scan statistic S m : 

P (S m (Lm) < n) « 201^02 (3 g) 

[l + qi - q2 + 2( qi -q 2 )^ L - 1) 

with an error bound of about 

E = {3 + r(a)(l - qi) + (L - 1) [1 + K(a)(l - g x )]} (1 - 9l ) 2 , (3.7) 
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where a = P(Wi > n) and qi,q 2 are defined by Eq.(1.2) as 

qi = P (Wi < n) = ¥(S m (2m) < n), (3.8) 

q 2 = P (Wi <n,W 2 <n) = P(5 m (3m) < n). (3.9) 

We should mention that in Haiman [1999, Theorem 4], the author obtained the 
following formula for the approximation error 

EH = {9 + 561(1 - qi ) + 3.3(L - 1) [l + 4.7{L - 1)(1 - qi ) 2 ] } (1 - q x f . (3.10) 

Notice that if qi, q 2 , q^, q^ are known then we can apply Theorem 2.5 to get a better 
approximation for the distribution of S mi even though in most applications Eq.(3.6) 
will suffice. 

Next we will restrict ourselves to the case of Bernoulli — 1 process, that is Y^s are 
Bernoulli trials with P(Yi = 1) = p = 1 — P(Yi = 0). In this particular framework, 
Naus (see for example Glaz, Naus and Wallenstein [2001], Theorem 13.1) provided 
exact formulas for P(S m (N) < n) when N = 2m and N = 3m, namely for q\ and 
<72 from Eq.(3.8) and Eq.(3.9). The following tables illustrates a compared study 
between the error formula given by Eq.(3.7) and the corresponding error bound used 
by Haiman [2007] along with the exact and approximated value for the distribution 
of the scan statistic: 



» 


qi 


q-2 


Approx 


Exact 


EH 


E 




Eq.(3.8) 


Eq.(3.9) 


Eq.(3.6) 




Eq.(3.10) 


Eq.(3.7) 


2 


0.97131 


0.95181 


0.82715 


0.82582 




0.01712 


3 


0.99716 


0.99500 


0.98001 


0.98000 


0.00032 


0.00010 


4 


0.99982 


0.99967 


0.99865 


0.99865 


1 x 10~ 6 


3 x 10~ 7 


5 


0.99999 


0.99998 


0.99994 


0.99994 


2 x 10~ 9 


6 x 10~ 10 


G 


1. 


1. 


0.99999 


0.99999 


1 x 10~ 12 


4 x 10~ 13 


7 


1. 


1. 


1. 


1. 


3 x 10~ 16 


9 x 10~ 17 



Table 3. Distribution of the scan statistic P (S m (Lm) < n) for 
m = 9, p = 0.05, L = 10 



» 


qi 


qi 


Approx 


Exact 


EH 


E 




Eq.(3.8) 


Eq.(3.9) 


Eq.(3.6) 




Eq.(3.10) 


Eq.(3.7) 


1 


0.96860 


0.94910 


0.74617 


0.74353 




0.02927 


2 


0.99813 


0.99677 


0.98061 


0.98060 


0.00019 


0.00006 


3 


0.99993 


0.99987 


0.99922 


0.99922 


2 x 10~ 7 


8 x 10~ 8 


4 


0.99999 


0.99999 


0.99998 


0.99998 


1 x 10~ 10 


4 x 10~ n 


5 


1. 


1. 


1. 


1. 


4 x 10" 14 


1 x 10- 14 



Table 4. Distribution of the scan statistic P(5 m (Lm) < n) for 
m = 10, p = 0.0165, L = 15 



The exact values for the distribution of the scan statistics presented in the tables 3 
and 4 (column " Exact" ) are computed using the Markov chain embedding technique 
described in Fu [2001]. 
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4. Proofs of the results 



4.1. Proof of Theorem 2.1. The proof will follow closely that of Haiman [1999]. 
Using the stationarity and 1-depcndcncc of the sequence we have 

p n = F(Xi > x, X 2 > x, . . . , X n > x) < V(X 1 > x, X 3 > x, . . . , X n > x) 

= F(X 1 > x)P(X 3 > x, . . . , X n > x) 

= PlPn-2, 

which give the basic inequality 

Pn<p [ ^ ] - (4.1) 

To show that C(z) has a zero in the interval (1, 1 + Ipi) it is enough to prove that 
C(l) > and C(l + Ipi) < 0. It easy to see that C(l) > 0, since 

oo 

C(l) = 1 + J2 i-^Pk-i = (Pi -P2) + ( P3 -P4) + --->0 (4.2) 

fc=1 >0 >0 

For C(l + Ipi) we have: 

oo 

C(l + l Pl ) = l + ^(-l) fc p fe _ 1 (l + Zpi) fc 
fe=l 

oo 

= -l Pl + V (1 + lpi) 2k \p2k-l ~ P2fe(l + Ipi)] 
k^l " " ' 

<P2k-l<Pi 

OO 

< -^i+^[(l + ipi)V] fc . (4.3) 

k=l 

It is easy to see that if (1 + Ipi) 3 < I then P i(l + Ipi) 2 < i+l pi < 1, so series in 
Eq.(4.3) is convergent and C(l + lp\) < 0. From the definition of I and the relation 
1 < t2(a) < t2(a) + e < we obtain (i 2 = t 2 (a)) 



4+t 2 {/t 3 2 + £+ V(*2 + e) < 4+t 2 (t 2 +e) + (t 2 +eY 

1 1 1 _ 1 
3a 3a 3a a 



which imply that t 2 — y/t 3 , + e + ae < 0. Combining this last relation with the fact 
that t 2 is a root of the equation at 3 —t + 1 = 0, we have that (1 + Ipi) 3 < l- For 
showing that the zero is unique we will prove that C'(z) is strictly decreasing on 
the interval (1, 1 + Ipi), i.e. C'(z) < 0. Using Lagrange theorem on [1, 1 + a] we get 
for 9 e (0, 1) 

0'(l + a)=C'{l)+aC"(l + 9a) (4.4) 
We will approximate both C"(l) and C"(l + 9a) as follows: 

C"(l) = -1 + 2pi - 3p 2 + R (4.5) 

where 

oo oo 

R = J2 2k (P2k-i -P2k)-^P2k- (4.6) 

fc=2 fe=2 
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To approximate R notice that £>2fc-i — P2k > 0, which implies 

oo oo 

- J^Pfflfe < R < 2tt (P2*-l ~ Kfc) 

fc=2 

and using Eq.(4.1) we get 



k=2 



Pi 



I -Pi 



= -J2p k 1 <R<2p l J^ kpl- 1 = 2p\ 



1 



(\-PiY i-pi 



k=2 k=2 

so that 

\R\ < 2p\ 

For C"(z) we have: 

C"{z) = 2pi - 2 • 3p 2 2 + 3 • 4p 3 z 2 - 4 • 5p 4 z 3 + 5 • 6p 5 2 4 



(i-pi) 2 i-pi 



(4.7) 

(4.8) 
(4.9) 



Using z £ (1,1 + Zpi) and Eq.(4.1) we have 

oo 

C"(z) < p 1 ^(2fc + l)(2fc + 2)(z^T) 



2 k 



fc=0 



< 2pi 



l + 3pi(l + /pi) 2 



[l- Pl (l + Z Pl )2] 3 ' 
Since z < 1 + Zpi, Eq.(4.9) and Eq.(4.1) implies 

oo oo 

C"(z) > -l Pl J2( 2 k + l)(2k + 2)p 2k+2 z 2k -4^2(k + l)p2 k+ 2Z 



(4.10) 



2fc+l 



fe=0 



> ^p 2 £ (2k + 1)(2* + 2)(z^T) 2 ' £ - 4 Pl z >T (* + l)(piz 2 ) fc 



fc=0 



fe=0 



> -2lpl- 1 + 3p ^ + lp ^ 



4 Pl (l + l Pl ) 



which in relation with Eq.(4.10) shows that 

l + 3pi(l + ? P i) 2 4pi(l + Zpi) 



\C"{z)\ < 2lp{ 



(4.11) 



(4.12) 



Combining Eqs.(4.4), (4.5), (4.8) and (4.11) we can show that C'(z) < if the 
following inequality is true 



- 1 



2pi 



(1-Pi) 2 



2lp{ 



2 i + 3 Pl (i + ; Pl ) 2 



[1 - P1 {1 + 1 P1 ) 2 } 



< 0. 



(4.13) 



We observe that the expression on the left hand side of Eq.(4.13) is increasing in 
Pi and since for p\ = 0.1 we have I < 1.1535 we get 



2pi 



(I -Pi) 2 

which verifies that C'(z) < 0. 



2lp\ 



l + 3pi(l + foi) s 



21 3 



< 0.3 



(4.14) 
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Now we try to approximate the zero A. From Lagrange theorem applied on the 
interval [1, A] we have C*(A) - C(l) = (A - l)C'{u), with u € (1, A) C (1, 1 + Ipi). 
Since C(A) = we get 

C « (4.15) 



A - 1 



C'{u) 



and taking /i = p\ — p 2 + P3 — pi + 2p 2 + 3p 2 ~~ 5pip 2 as in Haiman [1999] we obtain 

C(l) + pC'(u) 



A — (1 + p) 



C"(«) 



(4.16) 



Applying Lagrange theorem one more time as in Eq.(4.4) and combining with 
Eq.(4.5), the relation in Eq.(4.16) becomes 

C(l) + p(-l + 2pi - 3p 2 ) + p(-R + oC"(l + 0a)) 



A-(1 + m) = - 



C*(u) 



(4.17) 



where a = u — 1 and e (0, 1). 

If we denote T x = C(l) + p(-l + 2p a - 3p 2 ), then 

7l = (pi -p 2 +P3 -P4 +P5 -P6 + ...) _ M + ( 2 Pl - 3p 2 )p 

= (P5 - P6 + •••) + (Pi - Pa)(2pi - 3p 2 ) 2 + 2(pi - p 2 )(P3 - Pa) 
-P2{P3 - Pa- 
using Eq.(4.1) in Eq.(4.18) we observe that 

oo 

-P? < -P 2 (p 3 - P 4 ) < T! < ^ p* +1 + 4p? + 2p? 



(4.18) 



which gives us 



|C(l)+p(-l + 2pi-3p 2 )|<p? 



1-Pi 



(4.19) 



(4.20) 



Also we can notice that p > and 

P = (1 -pa)(Pi -Pa) +P3 -P4 + 2(pi -p 2 ) 2 
< (1 -p 2 )(pi -p 2 ) +Pi(pi -p 2 ) + 2(pi -p 2 ) 2 
= 3(pi -P2) 2 +pi -p 2 <Pi(l + 3pi). 
The last step is to find an upper bound for |C"(m)| _1 . For this observe that 
IC'iu)]" 1 = \l-(2p 1 u-3p 2 u 2 +4p 3 u 3 ~5p 4 u i + ...)\- 1 

< |l-|2piM-3p 2 u 2 +4p 3 M 3 -5p 4 it 4 + ...|p 1 (4.21) 

where we used the inequality |1 — x\ > |1 — \x\\. Taking into account that u £ 
(1, A) C (1,1 + ipi) and denoting the expression inside the second absolute value in 
the denominator of Eq.(4.21) by T 2 , we have 



To 



2k- 1 -P2kU)U 



2 k 



1 -J2 p2k ' 



,2k 



fc=l 



fc=l 



< 2{ Pl - P2 )uY,Kpm 2 ) k - 1 < Tl 



2p\u 



k=l 



(1 — PiU 2 ) 2 



(4.22) 
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In the same way 



T 2 > -^P2ku 2k - 2l Pl k P 2kU 2k ~ 1 



> 



k=l k=l 



k=l k=l 
.21 



T »\ < 7, \J. Lz - (4-24) 



2l Pl +u{l - Pl u 2 ) 2 Pl u 
> -p\u 5 > - — ^-r. (4.23) 

Combining Eq.(4.22) with Eq.(4.23) along with u < 1 + Ipi, we have 

2 Pl {l + l Pl ) 
[1 - Pl (l + l Pl ) 2 } 
Substituting Eq.(4.24) in Eq.(4.21) we obtain the bound 

I^MI -1 * ] 2M 1+IP1) ■ (4-25) 
[i- Pl (i+i Pl y]* 

Combining the Eqs.(4.8), (4.12), (4.19), (4.20), (4.25) along with the fact that 
\a\ < Ipi in Eq.(4.17), we obtain 

,,,,,, < |C(1) + //(-I + ggi - 3pg)j + + \a\\C"(l + 

|A - (1 + /i)l - \c^\ 

< K( Pl )p\ (4.26) 

where 

K ( pi ) = 1 2 Pl[ i+i Pl) ~ ■ (4-27) 

[l- Pl (l+Z Pl )^ 

To obtain K(a) it is enough to substitute p\ in the above relation with a with the 
additional remark that I = 1(a). □ 

4.2. Proof of Theorem 2.3. For completeness we will give a detailed proof of 
the theorem even if we repeat most of its ideas from Haiman [1999]. Remembering 
that 

oo oo 

C(z) = ckz k = l + J2 (-l) k Pk-iZ k (4-28) 

fe=0 fc=l 

we define 

oo 

D{z) = —- = Y, d ^ k (4.29) 
U[z) k=0 

which exists since c = 1, and from C(z)D(z) = 1 we have that do = 1 and 

n n 

djCn-j = (-lT~ j Pn-j-idj = 0, n > 1. (4.30) 

J=0 3=0 

Now if we define = {Xi <£,..., Xf. < x, > x, . . . , X n > x}, for k = 0, n, 

we obtain P(A ) = p„, P(A n ) = q„ and 

P(i4 fc )+P(^fc-i)=l^-fc9fc-i. fc >! (4-31) 
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Summing Eq.(4.31) over k, we deduce that 

n 

Qn = ^2 (-l)"~ fe Pn-fcft-i, n > 0, q_ 1 = q = l 



(4.32) 



k=0 



and comparing with Eq.(4.30), after using mathematical induction, wc conclude 
that d n +i = q n and that 



D(z) = ^2q k - 1 z k , g_i=gb = L (4-33) 

fe=0 

Taking A as in Theorem 2.1, we can write C(^) = U(z) (l — f) and observe that if 

n 

wc let U (z) = UkZ k then 



fc=0 



which shows that 



£ Ufc zM(i-f)=i + ]T(-i)W 



fc=l 



A 



= (-l)"p„_l, W = 1, ?1 > 1. 



Multiplying Eq.(4.35) with A™ and summing over n wc find 

n 

i + ]T(-i)V-iA fc 



fe=i 



A" 



-, n>l. 



If we denote with jj^y = 5Jifc2 fc = then 



fc=0 



D(z) 1 



A 



T(z) 



and using the same argument as above we get to = e?o = 1 and 

*n = «n — , n > 1 

so d n A™ = t + tiX + • • • + t n \ n , that is 

q n X n+1 =t +h\ + --- + t n X n + t n+1 X n+1 . 
To obtain the desired result wc begin by giving an approximation of u n : 



fe=i 



C*(A)=0 



53 (-l) fc p fc -iA A 

fc=n+l 



A" 



A'' 



< 



A" +1 



A" 



|i?n - Pn+lA + p n +2X 2 - Pn+3>? + ■■■ 
< X(\p n - p n +lX\ + \p n+ 2 ~ Pn+3A|A 2 + . . . ) 

Since A G (1,1 + Ipi) C (1, 1 + la) we have 

Pn - Pn+l(l + la ) <Pn~ Pn+lX < p„ - p n +l 

which shows that 

\Pn -p„+lA| < p n -p n+ l(l - la). 



(4.34) 
(4.35) 

(4.36) 

(4.37) 

(4.38) 
(4.39) 



(4.40) 

(4.41) 
(4.42) 
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If we denote by h = 1 — la and we use the bound from Eq. (4.42) and the fact that 
(Pn)n is decreasing in Eq.(4.40) we obtain 
\ u \ 

— ^- < Pn -Pn+lh + {Pn+2 ~ Pn+?,h)\ 2 + . . . 

= (Pn + Pn+2^ 2 + P«+4A 4 + ...)- h(p n+ i + p n+3 A 2 + p n+5 A 4 + . . . ) 

< W - h(p n+2 + Pn+A 2 + Pn+6^ 4 + ■■■) 

= W-^(W-p n ) = w(l-^)+^p n , (4.43) 
where by Eq.(4.1) 

W = p n + p n+2 >? + Pn+4^ 4 + ■■■ 

< p [ r 1] +p [ r 1] p 1 \ 2 + P [ r 1] pjx 4 +... 

= (1 + Pl A 2 +p 2 A 4 + ...)= x2 - (4-44) 

From Eq. (4.43) and Eq. (4.44) we conclude that 



'H ^ I 



1 *(l- 1 



1 -piA 2 A 2 V 1 - Pi A 2 



(4.45) 



l-PiA 2J 

Until now we have an approximation for u n , but we still need one for t n and to 
solve this aspect lets write 

^) = ^ = T3,rb(5))=g(-w-ir. (««> 

which is true since the convergence of C{z) implies \z\ < -^== so that 

\(1- Pl h) sr t| . n ^\z\X^pl[l- Pl (l-l Pl )} 



I—I <- ^X>?i^ d-^a-wvB) 



^T(l + ^i)(l + ^i) 2 f44?1 



Since uq = 1 we have U — 1 = M„z n and 



n>l 

(<7-l) fc = ]T (4-48) 

ij >lj — 

Combining Eq.(4.46) with Eq.(4.48) we get 

oo oo 

^^« = ^(-l) fe ^6^, (4.49) 

n>0 fe=0 (=1 

where 

k.k= ^2 u tl ...u tk . (4.50) 

iiH Mfc=_i 
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From Eq.(4.49) we get that to = 1 and 

n 

f„=^(-l) fe Yj U n ...U lk , k>l. 



k=l 



i x -\ Hfc_=n 

ij>X,j=X,k 



Notice that from Eq.(4.45) we can write 

iiH Yik—n 

ij>l,j=l,k 



n + ii 



•Pi 



where S = A jj_ p ^^ ■ Since, by induction it is easy to verify that 



k + 1 




"4 + 1" 


> 


'*i + 


■•■ + 4 + 1" 




"n + 1" 


2 


H h 


2 


2 




2 



(4.51) 



(4.52) 



(4.53) 



and the number of terms in the sum of Eq.(4.52) is equal with the number of 

different positive integers solutions of the equation i\ H h 4 = n, which is given 

by , we deduce 



U-i/' 



\b n ,k\ < Pi 

Now from Eq.(4.51) and Eq.(4.54) we have 



[n±l] (n 1 
k - 1 



(4.54) 



P^E f^V^Pf^E ("''W 



fc=0 



fc=0 



2fc + 1/ 



so that 



l*n|<rf[ Br] [(l + «) B - 1 + (l-<0 B - 1 ] 



We see from Eq.(4.39) and Eq.(4.56) that the difference 



\q n X n -q 3 X s \ = 



ra+l 



E**a* E'- A ' 



n+l 



s=5 



s=5 



(4.55) 



(4.56) 



r 00 fa+ll 

^ 5 EM 2 ^(l + ^ + a-^A- 1 (4.57) 



If we denote by o\ = (1 + <j)A, U2 = (1 — 5) A and by V the bound in Eq.(4.57), it 
is not hard to see that 



v= 5]ft\<4(l + tri) , ^(1 + ^2) 



1 - p\o{ 1 - picrj 



(4.58) 



Recalling that /i = 1 — Ipi, A € (1, 1 + Ipi) and p\ < 0.1 we observe that a?, is 
bounded by 

Ipi[l + 2pi(l + ipi)] 



l- Pl (l+l Pl f 



1-Pl ' 



DISTRIBUTION OF EXTREMES OF 1 DEPENDENT R.V. SEQUENCES 13 

which gives [0*2 ] < 0.5 and ^l-pt^) ^ Substituting the last relation in 

Eq.(4.58) we can rewrite the bound in Eq.(4.57) as 

\q n \ n -q 3 \ 3 \ <E( Pl )p 3 (4.59) 

where, if we denote by 77 = 1 + Ipx, 

, > V 5 [l + (l-2p 1 )r,f[l+p 1 ( V -2)} [i + v +(l-3 Pl )r,2] 

[Pl> 2{\- Pl ^y[{i- Pl ^y- PlV m + 7 1 -2 Pl r 1 Y] ■ [ - DUj 

To obtain E(a) it is enough to make in Eq.(4.60) the substitutions: pi a, 
I — > 1(a) and r) — > 1 + la. The following lemma gives an approximation for q 3 X 3 . 

Lemma 4.1. If pi < a then 

q 3 X 3 = 1- P2 + 2 P3 - 3p 4 + pf + G P 2 2 - &P1P2 + 0{L{a)p\) (4.61) 

where 0(x) is a function such that \0(x)\ < \x\ and L(a) is an ex P ression de P ending 
on a and given in Eq.(4-77). 

From Eq.(4.57) , Eq.(4.60) and Lemma 4.1 we conclude that 

\q n X n - (1 - p 2 + 2p 3 - 3p 4 + pj + 6p 2 2 - 6pipa)| < T(a)pl , (4.62) 

with T(a) = L(a) + E(a) and this ends the proof of the theorem. □ 

4.2.1. Proof of Lemma J^.l. From Eq.(2.1) of Theorem 2.1 we can write 

A = 1 + pi - P2 + P3 - Pi + 2 P \ + 3p 2 2 - 5p lP2 + 0(K(a)p\) (4.63) 
and raising to the third power we get 

A 3 = (l + Ci+(2) 3 + 3(l + Ci+C2) 2 0(K(a)pl)+0(K 3 (a)a 6 pl) + 

+3(1 + Ci + (2)0(K 2 (a)a 3 pl), (4.64) 
where we have used the following notations 

Ci = P1-P2, (4.65) 
C2 = p 3 -p4 + 2pl + 3p 2 2 -5pip 2 - (4.66) 

Since we can easily see from Eq.(4.65) and Eq.(4.66) that Ci = 0(p\) and (2 = 
0(3pl), we deduce that 

A 3 -(I + C1+C2) 3 = 0(S(a)p 3 ) (4.67) 
with S(a) given bellow by 

5(a) = 3(1 + a + 3a 2 ) 2 K(a) + 3a 3 (l + a + 3a 2 )K 2 (a) + a 6 K 3 (a). (4.68) 
Now we observe that by expanding (1 + Ci + C2) 3 we have 

(I + C1+C2) 3 = 1 + 3(Ci+C2 + Ci 2 ) + 6CiC2 + 3C 2 2 + 3CiC 2 2 + 3Ci 2 C2 

+C1 3 + cl 

= 1 + 3(Ci + C2 + Ci 2 ) + 0(18pl) + 0(27ap\) + 0(27a 2 p 3 ) 
+0(9a P 3 1 ) + 0( P 3 1 ) + (27 a 3 pi) (4.69) 
which along with Eq.(4.67) and Eq.(4.68) gives the following expression for 

A 3 = 1 + 3(Ci + (2 + Ci 2 ) + 0(P(a)pl), (4.70) 
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where P(a) can be computed by the formula 

P(a) = 3K(a)(l + a + 3a 2 )[l + a + 3a 2 + K(a)a 3 ]+a e K 3 (a) 

+9a(4 + 3a + 3a 2 ) + 19. (4.71) 
It is easy to see from Eq.(4.32) that 

q 3 = 1 - pi - 2(pi - pa) + pf - P3 = 0{{l - Pl ) 2 ) (4.72) 
and combined with Eq.(4.70) gives 

q 3 X 3 = q 3 [l + 3(Ci + C2 + Ci )] + 0(P(a)(l - Pl ) 2 p 3 ). (4.73) 

The last step in our proof is to find an approximation for the first term on the right 
in Eq.(4.73). If we write 

9.3 [1 + 3(Ci + C2 + Ci 2 )] = 1 - P2 + 2p 3 - 3p 4 + Pi + 6pl - 6piP2 + H, (4.74) 
then H checks the relations 

H = 3(p 1 ~p 2 ){(pl-p 3 )[3(p 1 -p 2 ) + l-p 2 ]+p 2 2 -9(p 1 -p 2 ) 2 } 
-3(P3 -P*)\pi +2(pi ~p 2 ) - (pj -p 3 )] 
= 0(36pl). (4.75) 

Finally, combining Eq.(4.73), Eq.(4.74) and Eq.(4.75) we get 

q 3 X 3 = l-p 2 + 2p 3 - 3p 4 + Pi + Gp 2 2 - &P1P2 + 0{L{a)p\) (4.76) 
where we used the notation 

L{a) = 36+ (1 - Pl ) 2 P(a) 

< 3K(a)(l + a + 3a 2 ) [1 + a + 3a 2 + K(a)a 3 } + a 6 K 3 (a) 

+9a(4 + 3a + 3a 2 ) +55. □ (4.77) 

4.3. Proof of Corollary 2.2 and Corollary 2.4. To prove the relation in Corol- 
lary 2.2 we see that 

1 +pi -p 2 + 2(pi -p 2 ) 2 =li- [P3 ~P±-Pi{P\ -Pi)\ (4.78) 

and since 

P3-P4 = P(Xi > x, X 2 > x, X 3 > x, X 4 < x) 

< p 1 ¥(X 1 >x,X 2 <x)=pl(p 1 -p 2 ), (4.79) 

we get 

|P3 ~P4-P2(Pl -P2)| < P3 - Pi +P2(Pl -Pi) < Pv (4.80) 
Combining these relations with Eq.(2.1) from Theorem 2.1. we obtain 

|A - (1 + Pl - p 2 + 2(pi - p 2 f)\ <p 2 (l + A»pi) < {l + aK{a))p 2 D (4.81) 

To prove Corollary 2.4 notice that 

- 3pj < P 2 + 2 P3 - 3 P4 + 6p 2 - 6 PlP2 < 3 P 2 (4.82) 

and using Eq.(2.4) from Theorem 2.3, we get 

\q n X n - (1 - p 2 )\ < T(a)p 3 + \p 2 1 +2p 3 -3p i + 6pl-6p lP2 \ 

< aT{a)p\ + 3p\ = (3 + aT{a))p\. □ (4.83) 
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4.4. Proof of Theorem 2.5 and Theorem 2.6. Denoting by 

Ml = 1 - P2 + 2p 3 - 3p 4 + p\ + 6P2 - 6piP2, 



/X 2 = 1 + Pi - P2 + P3 - Pi + 2Pl + 3^2 - 5 PlP2 



we observe that 
and that 



0<^<1 

M2 



fJ-2 = 1 + (pi -Pa)(l -P2) +P3 ~Pi + 2(pi -P2) 2 > 1. 



(4.84) 
(4.85) 



With the help of Eq.(2.1) and Eq.(2.4) we get 



Mi 
MS 



< 



A" 



A" A i» 



< r(o)p? + |Mi| 



1 

M2 




M 2 



< [r(a) + nX(a)]pf. 
To express [i± and //2 in terms of q's we have to observe first that 

Pi = 1-qi 
p 2 = 1 - 2qi + q 2 
P3 = l — 3qi + 2q 2 + q\ - qs 
Pi = 1 - 4gi + 3(? 2 - 2q 1 q 2 + 3g 2 - 2q 3 + q 4 
and after the proper substitutions, we get 

Mi = 1 + «i - ?2 + ?3 - qi + 2q\ + 3g<2 - 5gig 2 
M2 = 6(gi - q 2 ) 2 + 4gr 3 - 3<? 4 

To finish the proof of Theorem 2.5 it is enough to use the above relations in Eq. 
to obtain 

6(gi -g 2 ) 2 + 4g 3 -3 94 



(4.86) 



9» 



< Ai(l-gi) ; 



(1 + 9i - 92 + q 3 -Qi + 2q\ + 3g| - 5giga)' 
where 

Ai = Ai(a,n) = r(a) + nK(a).n 
For the proof of Theorem 2.6 we use the same approach as in Theorem 2.5. 
the help of Eq.(2.3) and Eq.(2.7) we get 



(4.87) 
(4.88) 
(4.86) 

(4.89) 

(4.90) 
With 



< 



A" 



vi_ 

X" 



i'i 



< (3 + piT(a))p 1 + nv\ 



|A-^ 2 | 
A v 2 



< [3 + T(a) Pl +n(l+piK{a))]pl 

where v\ = l—p 2 and v 2 = 1 +p\ — p 2 + 2(pi — p 2 ) 2 ■ We express v\ and v 2 in 
of g's using the relations for pi to p 4 from the proof of Theorem 2.5 and we 

vi = 2qi - q 2 

v-i = 1 + qi - 92 + 2(91 - q 2 ) 2 , 



(4.91) 

terms 
have 

(4.92) 
(4.93) 



16 



AMARIOAREI ALEXANDRU 



which wc substitute in Eq.(4.91) to obtain 




(4.94) 



A 2 = A a (a, n, ft) = 3 + r(a)(l - ft) + n [1 + #(a)(l - ft)] .□ 



(4.95) 
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